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Abstract 

The rare decay (fxj) can occur only via penguin annihilation topology in 

the standard model. We calculate this channel in the perturbative QCD approach. 
The predicted branching ratio is very small at (10^^). We also give the polarization 
fractions, which shows that the transverse polarization contribution is comparable to 
the longitudinal one, due to a big transverse contribution from factorizable diagrams. 
The small branching ratio in SM, makes it sensitive to any new physics contributions. 

1 Introduction 

The study of B meson decays has offered a good place to test the standard model (SM) 
and to give some important constraints on the SM parameters. Recently, more attentions 
have been paid to the B VV decay modes. The transverse polarization of the vector 
meson can contribute to the decay width, and the fraction of each kind of polarization has 
been or will be measured. In some penguin dominated decay modes, such as 5 — > 0i^r* 
PP, the experimental results for polarization differ from most theoretical predictions [2], 
which has been considered as a puzzle and lots of discussions have been given jHlH!- So the 
* Mailing address. 
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polarization problem in the B — > VV decay modes brings a new challenge to the standard 
model, maybe it is a signal of new physics j3J |3] . 

In this work we will calculate the branching ratio and the polarization fractions of the 
charmless decay channel 00 with perturbative QCD approach (PQCD) jnilZ]. In this 

channel, the initial b quark and the light valence d quark in the B meson don't appear in 
the final states, so it must be an annihilation topology in Feynman diagrams. Annihilation 
diagrams can't be calculated in factorization approach |H1 El or in QCD improved factoriza- 
tion approach ilOj for its endpoint singularity, but in PQCD approach this singularity can 
be regulated by Sudakov form factor and threshold resummation, so the PQCD calculations 
can give converging results and have prediction power. In this channel, since no tree level 
operators can contribute, the dominant contribution comes from penguin operators. The 
annihilation topology is usually suppressed relative to the emission topology which can ap- 
pear in other modes, so this channel is a rare decay mode, and hasn't been measured in the 
B factories. 

In the next section we give our theoretical formulae based on the PQCD framework. 
Then we show the numerical results and a brief conclusion in the third section. 



2 Perturbative calculation 

For simplicity, we work in the B meson rest frame, and adopt the light-cone coordinate 
system. Then the four-momentum of the B meson and the two mesons in the final state 
can be written as: 

Pi = ^(1,1,0t), 
P2 = ^(l-r,r,0T), 

P3 = ^(r,l-r,0T), (1) 



in which r is defined by r = i(l - y 1 - 4M|/M|) ~ M^/M| < 1. To extract the hehcity 
amplitudes, we should parameterize the polarization vectors. The longitudinal polarization 
vector must satisfy the orthogonality and normalization: " -P2 = 0, £3^^ ■ P3 = 0, and 
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= e^L = —1- Then we can give the manifest form as follows: 

e2L = ^=(1 -r, -r, Ot), 



e,L = ^{-r,l-r,OT). (2) 
\/2r 



As to the transverse polarization vectors, we can choose the simple form: 

e2T=^(0,0,lT), 

e3T = ^(0,0,lT). (3) 

Only penguin operators can contribute to this decay channel, so the relevant effective 
weak Hamiltonian can be written as [TTj : 

n,jf = ^Vt,v;,aif^)o,ifi), ^ = 3-10, (4) 

where Ci are QCD corrected Wilson coefficients, and Oi are the usual penguin operators 
with the form 

O3 = {Sibi)v-AY.iQjQj)v-A, Oa = {Sibj)v-AY.iQjQi)v-A 

Q q 

O5 = {Sibi)v-Aj2iQj(lj)v+A, Oq = {Sibj)v^Aj2iQj(li)v-A, 

q q 



O7 = l{Sibi)v~A E (^q{(lj(lj)v+A, Os = l{Sibj)v-A E (^q{<lj<li)v+A, 



(5) 



O9 = i^{sibi)v-A E (^q{(lj(lj)v-A, Olo = ^{sibj)v-A E eq{qjqi)v-A- 
q q 

where q = s. The first four operators are QCD penguin operators; while the last four are 

electroweak penguin operators, which should be suppressed by the coupling a/as- 

The decay width for this channel is: 

« (T=L,T 

where Pc is the 3-momentum of the final state meson, with |Pc| = ^(1 - 2r). Note that 
for our case an additional factor 1/2 should appear for the permutation symmetry of the 
identical final state particles. The decay amplitude A^'^ which is decided by QCD dynamics 
will be calculated later in PQCD approach. The subscript a denotes the helicity states of the 
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two vector mesons with L(T) standing for the longitudinal (transverse) components. After 
analyzing the Lorentz structure, the amplitude can be decomposed into |T]: 

= MIMl + M^MNelia = T) ■ el{a = T) + iMre.upA* € P2P^i ■ (7) 

We can define the longitudinal i^o, transverse H± helicity amplitudes as 

= MIMl. = MIMn T MIV?^Mt, (8) 

where r' = {P2 ■ P^) / M'^. After the helicity summation, we can deduce that they satisfy the 
relation 

^ W^M" = \Hq\^ + \H+\^ + \H.\^. (9) 

(T=L,R 

There is another equivalent set of definition of helicity amplitudes 

A\\ =^V2MIMn, 

= iMl^r'^ -IMt, (10) 

with ^ the normalization factor to satisfy 

|AoP+|A|||2 + |A^|2 = l, (11) 

where the notations Aq, A\\, A± denote the longitudinal, parallel, and perpendicular polar- 
ization amplitude. 

What is followed is to calculate the matrix elements A4l, and A4t of various op- 
erators in the weak Hamiltonian with PQCD approach. In PQCD approach, the decay 
amplitude is factorized into the convolution of the mesons' light-cone wave functions, the 
hard scattering kernel and the Wilson coefficients, which stands for the soft, hard and harder 
dynamics respectively. The transverse momentum was introduced so that the endpoint singu- 
larity which will break the collinear factorization is regulated and the large double logarithm 
term appears after the integration on the transverse momentum, which is then resummed 
into the Sudakov form factor. The formalism can be written as: 

J dxidx2dx3bidbib2db2b3db-iTr[C{t)^Bixi,bi)^^{x2,b2)^^ix3,b3) 
H{x,M,t)St{xi)e-^'^\ (12) 
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where the 6j is the conjugate space coordinate of the transverse momentum, which represents 
the transverse interval of the meson, t is the largest energy scale in hard function H, while 
the jet function St{xi) comes from the summation of the double logarithms In^Xj, called 
threshold resummation ^2]; which becomes large near the endpoint. 

The light cone wave functions of mesons are not calculable in principal in PQCD, but they 
are universal for all the decay channels. So that they can be constraint from the measured 
other decay channels, like B — > Kit and B —>■ tttt decays etc. [7J. For the heavy B meson, 
we have 

■i + MB)750B(a:,6). (13) 



For the longitudinal polarized meson 
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and for transverse polarized meson, 



[M^ ^2lM^)+ i2L n<Pl{x) + M^/0^(a;)], (14) 



In the following concepts, we omit the subscript of the meson for simplicity. 

Now the only thing left is the hard part H . In PQCD approach, it contains the corre- 
sponding four quark operator and the hard gluon connecting the quark pair from sea. They 
altogether make an effective six quark interaction. The hard part H is channel dependent, 
but it is perturbative calculable. When calculating the hard parts (shown in the Figure 1), 
the factorizable diagrams (a) and (b) have strong cancellation effects, which results in null 
longitudinal polarization contribution and null parallel polarization contribution. The per- 
pendicular polarization survives with a large factorizable contribution, which will be shown 
later to make a large transverse polarization. The detailed formulas with polarization M.Li 
M.Ni and M.T for each diagram are given in the appendix. According to PQCD power count- 
ing rules, the longitudinal nonfactorizable diagram should give the leading contribution, and 
the contributions from the other diagrams are suppressed by a factor r. 



5 
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Figure 1: Leading order Feynman diagrams for 



3 Numerical results and summary 



For the B meson wave function distribution amplitude in eq. pH|) . we employ the model 



Nbx'^{1 — xY exp 



xM 



B 



(16) 



where the shape parameter ujb = 0.4GeV has been constrained in other decay modes. The 
normalization constant Nb = 91.784GeV is related to the B decay constant J'b = 0.19GeV. 
It is one of the two leading twist B meson wave functions; the other one is power suppressed, 
so we omit its contribution in the leading power analysis TT. The (p meson distribution 
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amplitude up to twist-3 are given by 
3/^ 



J 6 



fT 

J 6 



x(l — x), 

3(1 - 2xf + 1.68C|(1 - 2x) + 0.69 



1 + (1 - 2x)ln 



X 



1 — X 



3/J 



3(1 - 2x)(4.5 - 11.2X+ 11.2x^) + 1.38 In 



X 



1 — X 



2V2iVe 

J 6 



xil — x) 



3/J 



4v^ 

with the Gegenbauer polynomials 



l + 0.2C|(l-2a;) 
^[1 + (1 - 2xf] + 0.24[3(1 - 2xf - 1] + 0.96C|(1 - 2x) 
1 - 22;)[1 + 0.93(10x2 - lOx + 1)], 



(17) 

,(18) 
(19) 
(20) 
(21) 
(22) 



(23) 
(24) 
(25) 

We employ the constants as follows [12]: the Fermi coupling constant Gp = 1.16639 x 
lO^^GeV^^ the CKM matrix element \VtlVtd\ = 0.0084, the meson masses Mb = 5.28GeV, 
= 1.02GeV, the decay constants = 0.237GeV, /J = 0.22GeV and the B meson 
lifetime r^o = 1.55ps. The results for the center value of the branching ratio is then 



and the helicity amplitudes are given by 



1.89 X 10" 



Ro = 0.65, Rn = 0.02, Rj 



0.33, 



(26) 



(27) 



which shows that the transverse polarization contribution is comparable to the longitudinal 
one. The relative strong phases, 0|j = arg (yly/Ao), 4>± = arg {A±/Aq) are given by 



198.34° 



195.48°. 



(2J 



Now we consider the contribution from different operators. In the factorizable diagrams, 
A^L = TWat = 0, because of the cancellation between diagrams of Figure 1(a) and 1(b). For 



Air, the QCD penguin operators O3, O4, O5 and Oq, contribute at the same level. In the 
nonfactorizable diagrams, the operator Oq give the most important contributions. If we omit 
the contribution from the electroweak penguin operators, the variation of the contribution 
from nonfactorizable diagrams (Figure 1(c) and (d)) is small, while that of the factorizable 
diagrams (Figure 1(a) and (b)) is large. The reason is that the electroweak penguin operator 
Og, which has a large Wilson coefficient, only presents in the factorizable diagrams. The 
overall contribution of electroweak penguin at the branching ratio level is less than 30%. 
We also test the contribution without twist-3 wave functions. We find that if we keep only 
twist-2 wave functions the total branching ratio doesn't change much, but the contribution 
from the factorizable diagrams will vanish, and the transverse polarization contribution then 
becomes very small. So the twist-3 wave functions give very important corrections to the 
polarization fractions. 

There are many theoretical uncertainties in the calculation. The next to leading order 
corrections to the hard part is a very important kind of uncertainty for penguin dominant 
decays. To test it, we consider the hard scale at a range 

max(0.75MBL'„, ^ < ta < max(1.25MBi:>„, ^, h, (29) 

max(0.75MBL>6, ^) < ^6 < max(1.25MBA, 7^, 7^), (30) 
02 03 02 03 

max(0.75MBF, 0.75MbDc, ^, ^) < tc < max(1.25MBF, 1.25MbDc, ^, (31) 

Oi 63 h 63 

max(0.75MsF,0.75Ms|X|^:^,:^) < ta < max(1.25MBF, 1.25Ms|X|5 1 1), (32) 

Oi 03 bi 03 

and other parameters are fixed. Then we can obtain the value area of the branching ratio as 

Br{B' <t><t>) = (1.891^11) X 10"^, (33) 

which is sensitive to the change of so the next to leading order corrections will give 
important contribution. The ratios |Ao|^, and are also very sensitive to the 

variation of i, because that the nonfactorized contributions decrease as the increasing of t, but 
the factorizable diagrams, which gives the main contribution of the transverse polarization, 
increase. The variety area of |Aop is about 0.41 — 0.81 . 
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Another uncertainty is from the meson wave functions, which is governed by other mea- 
sured decays [7J. The variation of the parameters will also give the corrections, such as the 
parameter Ub in the B wave function, if we assume its value area is 0.32 — 0.48, we will give 
the branching ratio 

Br{B° ^ 00) = (1.89l[]:l) X 10"^ (34) 

The ratios Rq, R\\, R± is not very sensitive to the change of Uh, because it only gives an 
overall change of branching ratio, not to the individual polarization amplitudes. 

In this paper, we calculate the rare decay channel 00 in PQCD approach and give 

its branching ratio and polarization fractions in SM. This decay occur purely via annihila- 
tion topology, and only penguin operators can contribute. We predict that it has a very 
small branching ratio of 10~^. This is so small that it will be sensitive to the new physics, 
such as supersymmetry etc. |HlEj, which may give a larger branching ratio. The current 
experiments only give the upper limit: Br{B^ cfxp) < 1.5 x 10~^ ^Tj, so the more accurate 
experimental results are needed to test the theory. 
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A factorization formulas 

In the factorizable diagrams, due to the identical particles at the final states cancellation 
occurs between the two diagrams figure (a) and (b). Only the perpendicular polarization 
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f oo poo 

Jo 



part survives, 

M"^ = -167rCj./BM| I dx2dx3 [ 

Jo Jo 

+ r{x2)r{xs){l + X3) + 0"(X2)0'^(X3)(1 + X3) + 0"(a;2)0"(x3)(x3 - 1)] S^itf 

{t)h{x2,X3,b2,h), (35) 

/•l /»00 f'QO 

l&T^CpfBMl / (ix2da:3 / h2dh2 / 63d63ras(t)[0"(^2)0"(a:3)(^2) 
Jo Jo Jo 

+0^(X2)0"(X3)(2 - X2) + 0«(X2)0"(X3)(2 - X2) + 0"(x2)0"(x3)(x2)]5^(i)' 

{t)h'{x2X^MM), (36) 

where the /i functions come from the integral on the transverse momentum, their manifest 
forms is 



h{x2,X3,b2,h) 



fHl{MBFh2)St{,x^)[e{h - b2)Jo{b2MBDa)H^{hMBDa) 



+e(b2 - bs)Jo{bsMBDa)H^{b2MBDa)], 



h'{x2, xs, 62, bs) = CYyH'o{MBFb3)St{x2)[e{bs - b2)Jo{b2MBD,)Hl(b^MBD,) 
+e{b2 - b2)Jo{bsMBD,)H^{b2MBD,)], 
with the notation F and D stand for: 

F = ■^/[(l - X2){1 - r) + Xsr][x3{l - r) + (1 - X2r)] 
Da = + r (1 - X3)] [1 - r(l - X3)] 



(37) 



(38) 



Db = v^[l - X2 + rx2]{l - rx2). 
t is the hard scale, which is chosen as 

ta = max(MBL»„, I/62, 1/63), h = max(MsA, I/&2, l/^s)- 
The Sudakov form factor is written as 



S^(t) = exp 



-.(a;2P2+, 62) - s{{l - X2)P2^ 62) - 2 T ^7(«.(/^'))l , 

J 1/62 ^^ 



(39) 



(40) 



(41) 



with the quark anomalous dimension 7 = — a^/Tr and the s{Q,b), the so-called Sudakov 
factor, which comes from the resummation of the double logarithms, is given as 



s{Q,b) 



i/b 



^(27i,-l-log2)+C^logg^ "^^^'^ 



. 67 TT^ 10 2^ , 7E 
+ ^¥-y-27^^+3^°^°ST 



a^(/x') 



TT 



1 ^ 

log — 



(42) 
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The nonfactorizable amplitudes for diagrams (c) and (d) are written as 

on i\j2 p roo 

Ml = -^-^l[dx]J bidhhdh(pB{xi){[-l+X2 + r{2-4x2)] 

0(X2)0(X3) + r(l + X2- X3)0^(X2)0*(X3) + r(-l + X2 + X3)(l)\x2)f (x^) 
+r(l -X2- X3)0"(X2)0*(X3) + r(3 -X2+ X3)(l)'(x2)(l)'{x3)}as{t) 

[C, + Ce - Cs/2 - Cio/2]/i,(a;i, X2, x^, h, h)S{t), (43) 

Mt = -^-^jjdx]J^ hdb,hdh<l>B{^i){[xs-^^3r]<j>{x2)H^3) 

+r(l - X2 + a;3)(/)*(x2)(/>*(x3) - r(l - X2 - X3)(p\x2)(p'{x3) 

+r(l - 0:2 - 2;3)0'(a;2)0*(a;3) + r(-l + X2 - a;3)0'(x2)0'(a;3)}as(t) 

[C4 + Ce - C8/2 - Cio/2](i)/i;(xi, X2, xs, h, h)S{t), (44) 

= y [c^a;] y hdhhdhMxi)r[-'2<p''{x2)r{x3) 

+(j>'^{x2)(f>'^{xs){l + X2- X3) - 2r{x2)r{x3)]as{t) 

[C4 + Ce - C8/2 - Cio/2] (t)/i„(a;i, 0:2, X3, 61, 63)^(t), (45) 

M% = ^^'^J^'^^^J b^db^b3dh3ct>B{xi)rf{x2)f{x3){l-X2 

+X3)as{t)[C^ + Ce - C8/2 - Cio/2](t)/i;(xi, X2, X3, b,, b3)S{t), (46) 

M^T = ^^^^^^g /" [^a;] f hidb,b3dh3(t)B{Xi)r[2(j)%X2)(j)\x3) 







v/6 

+0^(X2)0^(X3)(1 -X2- X3) + 20'^(x2)0''(x3)]«.(i) 



[C4 - Ce + C8/2 - Cio/2](i)/i„(xi, X2, X3, 61, &3)-5(i), (47) 
= _ 647rCVM| r^^j r bM^db3(l)B{xi)r(f^{x2W{x3){l-X2 
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-a;3)a.(t)[C4 - + C8/2 - Cio/2](i)/i;(a;i, X2, X3, 61, 63)^(^)- (48) 



The h functions are defined as 

K{X^,X2,X3MM) = -[e{b3-b^)Ub^MBF)Hl{b3MBF)+e{b,-b3) 

Jo{b3MBF)H^{biMBF)] Ko{MBD,bi), (49) 

ITT 

h',{xi,X2,X3,bi,b3) = -[e{b3-b^)Jo{b^MBF)H^{b3MBF) + eibi-b3) 

Ub3MBF)Hl{b,MBF)\ x{ ^ ' (50) 

\ KoiVXbi), X>0, 

11 



with the notations 



[dx] = I dxi / dx2 / dxs (51) 



^0 JO 



Dc = a/1 - [x2 - xi + r(l - X2- X3)][l - X3 - r(l - X2 - xs)] (52) 
= [x2 + xi - 1 + r(l - X2 - a;3)][a;3 + r(l - X2 - xs)]. (53) 



And the hard scale t is 



= maxiMeF, MbD„ l/h, (54) 
U = max(MBF, Mb V^i, V&a). (55) 

The Sudakov form factor is S{t) = SB(t)S^{t), with 

Ssit) = s(xiP+, 60 + 2 r ^7(a,(/.2)). (56) 

Ji/bi 
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